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Abstract 

We investigate the possibility of measuring the thermal Casimir force and its gradient in the 
configuration of a plate and a microfabricated cylinder attached to a micromachined oscillator. The 
Lifshitz-type formulas in this configuration are derived using the proximity force approximation. 
The accuracy for the obtained expressions is determined from a comparison with exact results 
available in ideal metal case. Computations of the thermal correction to both the Casimir force 
and its gradient are performed in the framework of different theoretical approaches proposed in 
the literature. The correction to the Casimir force and its gradient due to lack of parallelism of 
the plate and cylinder is determined using the nonmultiplicative approach. The error introduced 
in the theory due to the finite length of the cylinder is estimated. We propose that both static and 
dynamic experiments measuring the thermal Casimir interaction between a cylinder and a plate 
using a micromachined oscillator can shed additional light on the thermal Casimir force problem. 
Specifically, it is shown that the static experiment is better adapted for the measurement of thermal 
effects. 

PACS numbers: 31.30.jh, 12.20.Ds, 12.20.Fv, 77.22. Ch 
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I. INTRODUCTION 



The Casimir effect which was predicted more than 60 years ago [l| has recently attracted 
much theoretical and experimental attention (see the monographs {^-[(j] for a overview of 
the subject). The Casimir force originates from quantum fluctuations of the electromagnetic 
field. Thus, it is of the same origin as the van der Waals force, but acts at larger separations 
between the interacting surfaces where relativistic retardation effects become significant. 
The unified theory of both the van der Waals and Casimir forces was developed by Lifshitz 
7J, |8| . It describes the free energy and force of the van der Waals and Casimir interaction 
using the frequency-dependent dielectric permittivity of interacting bodies. 

The Casimir force acting between two neutral surfaces becomes dominant at separations 
below 1 fim. Until the last decade there was a lack of experimental information about this 
quantum phenomenon, but this situation has changed in recent years due to rapid progress 
in nanotechnology. Following two landmark experiments , many measurements of the 

Casimir and Casimir-Polder (atom-wall) forces have been performed with metal, semicon- 
ductor and dielectric surfaces (a review of all experiments can be found in Ref. 111]). 

The experimental and theoretical investigation of the Casimir force at nonzero temper- 
ature produced an unexpected result. It turned out that if conduction electrons in metal 
plates with perfect crystal lattices are described by the dissipative Drude model, this results 



in a violation of the Nernst heat theorem in the Lifshitz theory The same takes 

place for dielectric-type semiconductor and dielectric test bodies if free charge carriers, which 
are unavoidably present at nonzero temperature, are described by the Drude-type contri- 



bution to the dielectric permittivity 



tempe 

flfl 



14j-|l6| . The parameter range, where a violation 



occurs, was recently rederived and extended also to the geometry of a sphere above a plate 



171 ]. It was demonstrated that for metals with impurities the Nernst heat theorem in the 



Lifshitz theory is preserved 



1 8M2 lj] . However, in the case of metals with perfect crystal 



structure (which is the basic model used in condensed matter physics), and for dielectrics 



a satisfactory solution was not found. This even resulted in the attempts 22 



23 



to modify 
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3l|- 



the Lifshitz theory, and this led to a controversial discussion in the literature 
From the experimental side, it was shown that data exclude the Lifshitz theory, if the role 
of conduction electrons is taken into account by means of the d issip ative Drude model. This 



was repeatedly demonstrated for metal |3J2|436[], semiconductor [37J, |38| and dielectric [391.140] 
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test bodies. 

All modern experiments measuring the Casimir force (with the single exception of the ex- 
periment of Ref. 41] which was burdened by rather large error of about 15%) were performed 
using the configuration of a sphere above a plate. Bearing in mind that contradictions be- 
tween Lifshitz theory and thermodynamics, and experimental data discussed above, are of 
great concern for the foundations of quantum statistical physics 42J, it is worthwhile to con- 



sider the possibilities of alternative experiments. One more configuration of great promise 



that was discussed in the literature 
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44j is a cylinder parallel to a plane plate. In some 



sense it is intermediate between the configurations of two parallel plates and a sphere above 
a plate, by preserving some advantages of the latter while making the problem of preserving 
the parallelism less difficult than for two plates. An additional advantage is that the Casimir 
'orce in the cylinder-plate geometry was recently calculated exactly in an ideal metal case 



45 



471 ]. This allows one to reliably estimate possible errors introduced from the use of the 



proximity force approximation. 

The original proposal 44j considered a relatively large cylinder with length L = 2 cm 
and diameter 2R = 6.35 mm separated from the plate with a wide gap a > 1 /im. It has 



been shown 



491 ]. however, that for metal coated spherical lenses with centimeter-size 



curvature radii, the electrostatic calibration of Casimir apparatus meets serious difficulties. 
Specifically, the residual potential difference may become dependent on separation, and the 
force- distance relation for the electric force deviates from the form predicted by classical 
electrodynamics in a sphere-plane geometry. (This does not happen in the experiments 
of Refs. 32l-l36| dealing with spheres of about 150 fim in radius interacting with the plate 
of an oscillator.) Similar anomalies were reported in Ref. 501 for the configuration of a 

n 

cylindrical lens above a plate with lens parameters L = 4 mm, R = 12 mm. In Ref. [51[ 
this anomalous behavior of the electric force was attributed to unavoidable deviations of 
the mechanically polished and ground surfaces from the perfect spherical shape assumed in 
elecrodynamical calculations. As recognized in Ref. js(3], the same effect might be responsible 
for the calibration problems arising in the case of a centimeter-size cylinder above a plate. 
Because of this, the consideration of much smaller cylinders seems to be preferable. 

In this paper we investigate the possibility of combining the advantage of high precision 
provided by the setup of a micromachined oscillator and of a cylinder-plate configuration 
with a cylinder radius of about 100 /im. For this purpose we consider the thermal Casimir 
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interaction between a metallized plate and a microfabricated metal coated cylinder attached 
to a micromachined oscillator. The potentially increased precision of force measurements 
combined with the significantly decreased sizes of cylinders, as compared with Refs. 44, 5o| . 
requires computations of the thermal Casimir force with a more sophisticated account of 
such factors as the lack of parallelism of the cylinder and plate and the finiteness of the 
length of a cylinder. Hncre, we derive the Lifshitz-type formulas for the Casimir force and 
for its gradient for a cylinder and a plate made of real materials using the proximity force 
approximation (PFA). The accuracy of the results is determined by a comparison with the 
exact expressions available for a cylinder and a plate made of ideal metal [46] . Computations 
of the thermal correction to the Casimir force and to its gradient using the resulting formulas 
are performed in the framework of different theoretical approaches proposed in the literature. 
The correction to the Casimir force and its gradient due to the lack of parallelism of a plate 
and a cylinder is determined in the nonmultiplicative way. The error due to the finiteness of 
the length of a cylinder is estimated using the results obtained with the help of world-line 
numerics 52| . Both static and dynamic experiments measuring the Casimir force and its 
gradient in the cylinder-plate configuration using a micromachined oscillator are proposed. 
It is shown that the static experiment is better adapted for the measurement of the thermal 
correction to the Casimir force. 

The paper is organized as follows. In Sec. II we derive the Lifshitz-type formulas for a 
cylinder-plate configuration. Section III is devoted to the computations of thermal correc- 
tions to the Casimir force and to its gradient. The correction due to the lack of parallelism 
of a cylinder and a plate is found in Sec. IV. In Sec. V we estimate the error introduced by 
the finiteness of a cylinder. The proposed experiments are discussed in Sec. VI. Section VII 
contains our conclusions and discussion. 



II. LIFSHITZ-TYPE FORMULAS FOR A CYLINDER ABOVE A PLATE MADE 
OF REAL MATERIALS 

Let us consider the upper surface of a metallized plate nearest to the cylinder as a 
coordinate plane z = 0. For convenience of calculations, in Sees. II-V we consider the 
cylinder arranged above the plate. All the results obtained, however, refer equally well to 
the cylinder arranged below the plate (see Sec. VI). In real experiments, the cylinders used 
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may have elliptical cross-sections. Elliptical cylinders can be considered using the same 
methods, as applied below to circular cylinders. The cylinder axis (coinciding with y axis) 
is parallel to the plane z = at a separation distance R + a, where R is the cylinder radius, 
and a is the minimum separation between the cylinder surface and the plate. The lower half 
of the cylinder surface is described by the equation 



z(x) = R + a- VR 2 - x 2 . (1) 

We approximately replace this half with a set of thick plane plates extended along the y 
axis made of the same material as the cylinder and parallel to the oscillator plate. Let us 
suppose that the Casimir pressure P (z(x),T) between each of these plates and its projection 
on the plane z = at temperature T is known. According to the Derjaguin method 6|, |53| . 
under the condition a«i? the Casimir force acting between a cylinder and a plate can be 
approximately calculated as 

F{a,T)= P(z(x),T) dxdy, (2) 

J cr 

where integration holds over the projection of a cylinder on the plane z — 0. Incorporating 



some additional assumptions, the Derjaguin method was reformulated as the PFA 



54| and 



widely used in the literature for the calculation of the Casimir 



.0 



orce between a sphere and 



|55|). 



a plate (different versions of this method are discussed in Refs. 

We assume that the material of the plate and the cylinder (Au) is characterized by the 
frequency-dependent dielectric permittivity e(u). (In fact, as described in Sec. VI, both 
the plate and the cylinder are coated with an Au layer with sufficient thickness such that 
they can be considered as made of bulk Au.) The Casimir pressure between the plates at a 
separation z(x) at temperature T is given by the Lifshitz formula 



k t 00 / r°° 



1=0 

r 2 e -2q l z(x) 

X 1 ^ r 2 e -2q lZ (x)- ( 3 ) 

a a 

Here, ks is the Boltzmann constant, k± is the projection of the wave vector on the plane 
z = 0, q t = (k\ + £f7c 2 ) 1/2 , and f, = 2nk B Tl/h with I = 0, 1, 2, . . . are the Matsubara 
frequencies. The primed summation means that the term with / = is multiplied by 1/2. 
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The reflection coefficients r a for the two polarizations of the electromagnetic field (transverse 
magnetic with a = TM and transverse electric with a = TE) are given by 

£iqi - h 



r TM = r TM (z&, k A 



mi + h ' 

qi - h 



r T E = r TE (i£ h k±) = 1 / , (4) 

qi + h 

where ki = [k\ + e^f / c 2 ] 1 ^ 2 and Si = e(i^i). Notice that Eq. fl3]) can be identically repre- 
sented as 

k T 00 i f°° 

00 

J2J2 r - ne ~ 2nqiz(x) - ( 5 ) 



X 

n=l a 



Now we substitute Eq. fl5]) into Eq. ([2]) and arrive at the expression 

F(a, T) = - 2keTL y ! qi k ± dk ± 
z=o 1/0 



^ pli 

X y"( r TM + r TE) / to" 2 ™ 9 ^ 
n=l ^° 



(6) 



where L is the length of the cylinder which is assumed to be infinitely large. Using Eq. ([I]) 
we replace the integration with respect to x with an integration with respect to z. We then 
introduce two new integration variables, u = (z — a) / R and v = 2g ; a, and rewrite Eq. (jSj) in 
the form 



00 



1=0 ° n 

-1 q _ u ) e -n§vu 



pi 

1 ./0 



X 

n=l 



1- (1-u) 2 

Here, the dimensionless quantity r is defined as r = 47r fc^Ta/^c), and the reflection co- 
efficients are expressed in terms of dimensionless variables v and Q = £;/w c = rl, where 
w c = c/ (2a), in the following way: 

eiv - vV + (ij - i)Cf 



r TM = r TM (w,0) 



r TE _ r TE (v, - — — ======. (SJ 

v + 2 + (e, - 1)C/ 
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Within the PFA, only the term of leading order in a/R is physically meaningful. Because 
of this the integral with respect to u in Eq. © can be evaluated as 



Jo v 71 - 
1 

due- n * vu 



na 



+ 



2nvR 

Substituting Eq. fl9]) into Eq. ([7]) we obtain 



2u 

a\3/2 

r) 



+ 0(u 



3/2 N 



(9) 



F(a,T) = 

oo 

E 



k B TL R 



4v/?a 2 V 2a 



E' 

1=0 



?l 2 dv 



X 



n=l 



„2n 
TM 



+ r 



2n \ 
TEJ 



(10) 



Using the definition of the polylogarithm function, the Lifshitz-type formula ( FlUj) describing 
the Casimir force for a cylinder in close proximity to a plate can be represented in the form 



F(a,T) 



k B TL R 



A^Tta 2 V 2a 



E' 

1=0 



, z ' 2 dv 



X 



[Li 1/2 (r| M e- v ) +U 1/2 (r 2 TE e- v )] . 



(11) 



This equation can be used in numerical computations of the Casimir force acting between a 
cylinder and a plate (see the next section). 



In the case of a cylinder and a plate made of ideal metal, r. 
temperature Eq. (TlT!) takes the form 

k B TL I R 



2 

TM 



TE 



,-iIM, 



o,0) 



dl 



,3/2 



dvj2 



1. Then at zero 



(12) 



2V?ra 2 V 2a j j rl ^ y/n 

Introducing the new variable £ = tI, changing the order of integrations and calculating the 
integrals with respect to ( and v one obtains 

IhhcL IR ^ 1 n 3 hcL IR 



a,0) 



647ra 3 V 2a ^— ' n* 1 



384a 3 V 2a' 



(13) 



Equation ffTB"]) coincides with the familiar expression for the confi gur ation of an ideal metal 
cylinder above an ideal metal plate obtained using the PFA in Ref. (see also Refs. jsiQ)- 
Both Eqs. (11 ip and (fT3|) are approximate results. Bearing in mind that Eq. (Ilip can be 
used for the comparison of experimental data with theory, it is important to estimate how 
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accurate it is. Fortunately, for an ideal metal cylinder above an ideal metal plate at zero 
temperature the exact analytical result for the two leading terms of the expansion of the 
Casimir force in powers of a/R is available [46] 

*?M) = -^t/£(i-Cu£ 



384vra 3 V 2a V eA i? 

c ** = 4 - -L ~ °- 2886 - ( i4 ) 

7T Z 60 



This result was confirmed by means of numerical computations in Ref. [471 ] . From the 
comparison of Eqs. ffl3|) and ( I14p it can be seen that under the condition a <^ R the 
deviation of the PFA result from the exact result does not exceed 0.3a/ R. For a typical 
experimental value oia/R ~ 10~ 3 the error in using the PFA turns out to be of about 0.03%. 
For the configuration of an ideal metal sphere above an ideal metal plate at temperature 
ksT 3> he/ R it was shown 56] that in the zeroth order of a/R the PFA leads to exact 
results for the zero-temperature contribution to the Casimir force and, separately, for the 
thermal correction to it. It was also demonstrated that at zero temperature for a sphere and 
a plate made of real metals the deviation of the PFA results from the exact results is smaller 
than for a sphere and a plate made of an ideal metal 57H60I]. Finally, exact computations of 
the thermal Casimir force at nonzero temperature in a sphere-plate configuration performed 
with different a/R show that the exact results approach the PFA results when a/R — > 



60 



62| . All the above suggests that for the configuration of a real metal cylinder above a 



real metal plate, the error of the Lifshitz-type Eq. (ITTj) obtained using the PFA is of about 
the same, 0.3a/ R, as holds for ideal metal bodies at T = 0. 

The Lifshitz-type formula for the gradient of the Casimir force can be obtained in a 
similar way to Eq. (fffj) . For this purpose we differentiate Eq. (JHJ) with respect to a taking 
into account Eq. ([[]), 



oo 



71 1=0 



°° pR 

xJ2 n ( r ™ + r T n E) dxe- 2nqiz{x) . (15) 

n=l J 
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Dealing with this equation as with Eq. fl6]), in the leading order of a/R one obtains 

dF(a,T) k B TL [R ' 



da 4v / 7ra 3 , _ 

v 1=0 

oo 



-Y v^dv 



X)v^e- nw (r& + r&). (16) 



x 

n=l 



Using the definition of the polylogarithm function results in 

dF(a, T) k B TL Hr 



^ poo 

V / v 5/2 dv 
[LU /2 (4 M e- v ) + Li_ 1/2 (4 E e-")] . (17) 



For an ideal metal cylinder above an ideal metal plate at zero temperature Eq. (TTTj) leads 
to the result 

dF 1M {a,0) _ 77r 3 hcL FR_ 
da 768a 4 V 2a' 

which is also immediately obtainable by the differentiation of Eq. (I13j) with respect to a. 

In the high-temperature (large-separation) limit only the zero-frequency terms of 

Eqs. pip and (TTTj) contribute to the result. Specifically, for the ideal metal case one obtains 

from Eq. (PIT]) 



4^a 2 

3((3)k B TL IR 



p /-oo 

£ / ^/ 2 Li 1/2 (e-) <fc 
2a i 

(19) 



16a 2 V 2a' 

where £(2) is the Riemann zeta function. For the gradient of the Casimir force this results 
in 

dF lM (a,T) = 15C(3)fc B TL Hr 

da 32a 3 V 2a' [ ' 

For real metals the asymptotic behavior at high temperature depends on the assumed model 

of dielectric properties (see the next section). 

Although the immediately planned experiments are for metallic test bodies, the case of 
a dielectric cylinder above a dielectric plate is also of some experimental interest. For 
dielectrics with a finite dielectric permittivity at zero frequency, e(0) = £q, the high- 
temperature behavior of the Casimir force is simply obtainable from Eqs. f lTO"]) and fTTTj) . 
In this case at zero frequency we have 

r T M(0,v) = =r , r TE (0,t;)=0. (21) 
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The zero-frequency term of Eq. (flOj) then takes the form 

Calculating the integral with respect to v the high-temperature behavior of the Casimir 
force between a dielectric plate and dielectric cylinder is obtained: 

pdiel/ rr-T\ 3k B TL R ^r$ n 
F ( °' T) = "^2^V2^^^ 

n=l 



31cbTL I R T . 2 



T^LisW). (23) 



32a 2 V 2a 

The same equation, where r 2 , is replaced with r , describes the high-temperature Casimir 
force acting between a dielectric plate and metallic cylinder (or, equivalently, metallic plate 
and dielectric cylinder). The latter result does not depend on the model used for the 
description of dielectric properties of the metal. 



III. COMPUTATIONS OF THE THERMAL CORRECTION TO THE CASIMIR 
FORCE BETWEEN A CYLINDER AND A PLATE COATED WITH GOLD 

Here, we use Eqs. ( TTTj) and (ITT)) for computations of the relative thermal corrections to 
the Casimir force and to the gradient of the Casimir force defined as 

w _ F(a,T = 300 K)-F(a,0) 
6t ~ F(a,T = 300 K) ' [M) 
(2) = F'(a,T = 300 K)-F' '{a, 0) 
T ~ F'(a,T = 300 K) 

This requires knowledge of the complete data for the dielectric permittivity of Au along the 
imaginary frequency axis. Usually these are found using the dispersion relation 



2 rwlme(u) , 



where Im£(o;) is calculated from real and imaginary parts of the complex index of refraction 



n, tabulated, e.g., in Ref. 63]. Keeping in mind that Ren(w) and Imn(u) are known only 
within a restricted frequency range (from 0.125 eV to 10000 eV for Au in Ref. 63J) the 
problem of extrapolation of Ime(u) to lower and higher frequencies arises (in fact only the 
extrapolation to lower frequencies is of practical importance). According to the Drude model 
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approach [6l. I 111 ll8H2lj. the extrapolation of Ime(u) to lower frequencies is performed using 



the imaginary part of the dielectric permittivity 



e D (u) = 1 - /V (26) 

CU{U! + «7) 



where u p is the plasma frequency and 7 is the relaxation parameter (for Au we use 64] 
tu p = 9.0 eV and 7 = 0.035 eV). 

According to the plasma model approach ay , |36|, |65j, |66| , the dielectric permittivity 
can be expressed in the following oscillator form 

2 N 

sM = l--f + X; 2 f ■ , (27) 

where the parameters of oscillators Uj 7^ 0, jj, and gj are determined [3] from the best 
fit of lme(uj) to the tabulated optical data for n related to interband transitions of core 
electrons. As discussed in Sec. I, the Drude and the plasma model approaches lead to 
drastically different predictions for the thermal correction to the Casimir force in the sphere- 
plate geometry. Here, we perform all computations using both approaches in order to find 
additional evidence which could help to resolve this difference. Note that simple Drude 
and plasma models, which do not take interband transitions into account, were used in 
Ref. {43} to compute the Casimir interaction in a cylinder-plate configuration at nonzero 



temperature. Recently an alternative approach was also proposed 67]. It aimed to determine 
e(i^) from both Ime(u) and Kee(u) within a frequency range where they are measured 
without using any extrapolations, but assuming the Drude-type behavior at zero frequency. 



It was shown [68[, however, that large errors in the measurement data for Ree(w) complicate 
the application of this approach in practical computations. 

In Fig. 1 we present the computational results for the quantities S^ D (line 1) and 8^ D 
(line 2) as functions of separation in the region from 0.1 to 5 /im. The computations were 
performed using the Drude model approach. As can be seen in Fig. 1, the thermal corrections 
to both the Casimir force and its gradient in a cylinder-plate configuration are negative over a 
wide range of separations. This is typical for a Drude model approach in the configurations 
of two parallel plates and a sphere above a plate as well. It should also be noted that 
the relative thermal correction achieves rather large values at short separation distances. 
Specifically, at experimentally relevant separations a = 150, 200, 300, 500, and 750 nm the 
thermal correction to the Casimir force is equal to -1.8%, -2.7%, -4.6%, -8.6%, and -13.9%, 
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respectively. At the same respective separations the thermal correction to the gradient of 
the Casimir force achieves more moderate values -0.1%, -1.6%, -2.9%, -5.6%, and -9.3%. 
The magnitude of the relative thermal correction to the Casimir force \5^ D \ achieves its 
maximum value 41.6% at a = 2.55 /im, while the maximum magnitude of the thermal 
correction to the gradient of the Casimir force \5^ D \ = 52% occurs at a = 3.6/im. Thus, 
at large separation distances, dynamic experiments, where OF /da is a directly measured 
quantity, are most suitable for the detection of the thermal correction predicted within the 
Drude model approach. 

In the high-temperature (large-separation) limit it is easy to obtain the analytic expres- 
sion for the thermal Casimir force calculated using the Drude model approach. Taking into 
account that Eqs. (jSj) and (I2"6"j) lead to r TM (0,t>) = 1 and r TE (0,t>) = 0, the contribution of 
the zero-frequency term in Eq. ffTTj) is given by 



F D (a,T) = -|^r/f f >Li 1/3 (0 *> 
8V7ra 2 V 2a J 



32a 2 V 2a 

For the gradient of the Casimir force 

dF D (a, T) 15((3)k B TL [R 



(28) 



(29) 



da 64a 3 V 2a 

is obtained. 

We next present the computational results in the case when the plasma model approach is 
applied to describe metallic coatings on a cylinder and a plate. In Fig. 2(a) the computational 
results for the quantities 5^ (line 1) and 8^ p (line 2) are shown as functions of separation in 
the region from 0.1 to 5 /im. In Fig. 2(b) the same results in the region from 0.1 to 1 /xm are 
shown at a larger scale. As for the configurations of two parallel plates and a sphere above 
a plate, the thermal correction for a cylinder above a plate is positive when the plasma 
model approach is used. As can be seen in Fig. 2(a,b), at short separations the relative 
thermal correction obtained within the plasma model approach is negligibly small. Thus, 
the correction to the Casimir force is equal to 0.016%, 0.024%, 0.044%, 0.13%, and 0.38% 
at separations a = 150, 200, 300, 500, and 750 nm, respectively. At the same respective 
separations, the thermal correction to the gradient of the Casimir force is equal to 0.0014%, 
0.0023%, 0.0046%, 0.012%, and 0.029%, i.e., it is even less. This makes the observation of 
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the thermal correction at short separations, as predicted by the plasma model, impossible 
in either the static or dynamic experiments. However, at larger separations the thermal 
correction to the Casimir force computed using the plasma model approach quickly increases. 
It is equal to 0.9%, 7.2%, and 46% at separations a = 1, 2, and 5 /xm, respectively. At the 
same respective separations, the thermal correction to the gradient of the Casimir force is 
equal to 0.063%, 0.82%, and 26.7%. Thus, within the plasma model approach the relative 
thermal correction to the Casimir force is always larger than to the gradient of the Casimir 
force. 

As in the case of the Drude model approach, one can obtain analytic expressions for the 
thermal Casimir force and its gradient at high temperature. From Eqs. (jS]) and ( 12 7p . for the 
reflection coefficients at zero frequency, one obtains 

r TM {0,v) = l, (30) 



. OLD y kjl u i x , I — — — — T\9 

rTE{0,v) = = -{av - Va 2 v 2 + 1) , 

av + y a z v z + 1 

where a = oj c /uj p = 5 / (2a). The quantity 5 = c/u p has the meaning of the skin depth in 
the frequency region of infrared optics. Equation ( ITTjl then leads to 

F p (a,T) = -p^M Tv^dv (31) 



— y 1 a 2 v 2 + 1 



8^a 2 V 2a 
x [Li 1/2 (e~*)+Li 1/2 (4 E (0,v)e- v )]. 

The integral of the first term on the right-hand side of of Eq. (!3~T|) was calculated in Eq. (|28|) . 
The integral of the second term can be calculated approximately with the help of perturba- 
tion theory in powers of the small parameter a (for Au 5q ~ 22 nm a). Using Eq. ( 1301 ). 
one obtains 

r 2n (o, v) = l- Anva + 8n 2 v 2 a 2 + 0(a 3 ). (32) 
The second integral on the right-hand side of of Eq. ( |3T|) is then found using this equation: 

; 3 / 2 Li 1/2 (r 2 TE (0,v)e- v )dv 



oo 

V 



oo _ 

E4 



, Jn 

n=l 



oo 

v 3/2 (l - Anav + 8n 2 a 2 v 2 )e~ nv dv 



3AC(3) (1 _ 1Qa + 70a2 
3AC(3) 



5 35 f 5 x 1 
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By adding the contribution of the first integral in Eq. ( )3~T]) . we obtain the high-temperature 
asymptotic behavior of the Casimir force 



F p {a,T) 



3((3)k B TL jR 
V 2a 



16a 2 



^5^0^35/50 
2 a 4 la 



(34) 



The high-temperature asymptotic behavior for the gradient of the Casimir force is given by 



dF p (a,T) _ 15({3)k B TL jR_ 

V 2a 



da 



32a 3 



1 _ 7_5_o + 63 (So 
2 a 4 la 



(35) 



From the comparison of Eqs. fl34"j) and fl35l) with Eqs. ff28|) and ff29l) . respectively, it can 
be seen that under the condition a -C R the main contributions to the high-temperature 
Casimir force and its gradient calculated using the plasma model approach are two times 
larger than the Casimir force and its gradient calculated using the Drude model approach. 
The same holds in the configuration of two parallel plates and for a sphere above a plate 



under the condition b « ^ 111]. 



IV. CORRECTION DUE TO LACK OF PARALLELISM OF A PLATE AND A 
CYLINDER 



Experimentally it is impossible to preserve the exact parallelism of a plate and a cylinder 
as was always assumed in the above calculations. Because of this, it is important to find 
the correction to the derived results which arises if there is some small angle 9 <C 1 between 
the plate and the axis of a cylinder. For ideal metal bodies such a correction factor which 



depends on the mean separation between the cylinder and plate was found in Ref . [44[ . The 
magnitude of the Casimir force between real metal cylinder and a plate taking into account 
their nonparallelity was also estimated by multiplying the result obtained for parallel bodies 
by the same factor. Calculations show, however, that for real metal bodies corrections 
due to deviations from perfect geometry (for instance, due to surface roughness or due 



to a nonparallelity of two surfaces) are not of precisely multiplicative character |6l. 111]. 
Specifically, at short separations (typically at a < 200 nm) the multiplicative approach may 
lead to large errors. Because of this, we here calculate the thermal Casimir force between a 
cylinder and a plate in a more inclusive way by applying the PEA to account not only for a 
cylindrical geometry but for a nonparallelity of a cylinder and a plate as well. 
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We start with Eq. for the Casimir force between a plate and a cylinder parallel to 
it. It is necessary to take into account that the cylinder is now inclined by a small angle 9 
with respect to the plane z = 0. This means that thick plane plates extended along the y 
axis, which are used in Sec. II as a substitution for a cylinder in the PFA, are also inclined 
at an angle 9. We replace each of these inclined plates with a set of plates of area dxdy 
parallel to the plane z = 0. This is equivalent to replacing z(x) in Eq. ([6]) with z(x) + 9y, 
and subsequently averaging over the length of the cylinder. The Casimir force then takes 
the form 



where a is now the mean minimum separation between the cylinder and the plate. At 9 = 
Eq. fl36|) coincides with Eq. (Q. 

The integration with respect to y in Eq. (1361) is trivial. We then introduce the other 
integration variables u and v, defined in Sec. II, and perform calculations similar to those 
leading to Eq. ( ITUl) . Finally, retaining only the leading order contribution in the small 
parameter a/R, we arrive at 



where one additional small parameter is defined as Aq = 9L/(2a). For ideal metals at zero 
temperature Eq. (I36I) reads as 




L/2 



(36) 




(37) 




oo 





oo 



(38) 



After performing the integration and summation, Eq. (!38|) results in 



F IM (a,0,9) 



K(A e )F 1M (a,0), 



(39) 
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where F lM (a,0) holds for an ideal metal cylinder parallel to an ideal metal plate, and is 
defined in Eq. ffl3|) . and the correction factor due to nonparallelity is given by 

1 



k(A 6 



5A 



_(1-As)*/* {l + A e fl\ ' 



(40) 



Note that this correction factor was first obtained in Ref. (44j by direct application of the 
PFA to an ideal metal cylinder inclined at an angle 9 to an ideal metal plate. Here we have 
shown that Eq. ( flOl) follows from a more general Lifshitz-type formula ( 1381) which takes into 
account the nonmultiplicative effects due to nonparallelity of a plate and a cylinder. The 
correction factor n(Ag) can be also used at nonzero temperature in order to estimate the 
Casimir force between a cylinder and a plate under an angle 9 by means of the approximate 
multiplicative approach 

F m (a,T,9) = K (A e )F(a,T). (41) 

Here we estimate the accuracy of this approach at different angles and separation distances. 

The same procedure, as described above for the Casimir force, can be used in order to 
find a nonmultiplicative effect of a nonzero angle between a cylinder and a plate on the 
gradient of the Casimir force. As a result, Eq. ffl6|) is replaced by 



dF(a,T) k B TL f~R 



D °° , POO 

2a Jn 



da 4-i/vra 3 , - 

v 1=0 

xE^-(4 + 4)^. (42) 

n=l 

Using Eqs. (ITT)) and (157)) we performed numerical computations of the quantity 

^\a,A e ,T) = ^^l (43) 

at different a and Ag. This quantity takes into account the nonmultiplicative effects due to 
a nonzero angle between a cylinder and a plate. When compared with k(A#), it allows us to 
determine the application region of the multiplicative approach. Note that ft( nm ) depends on 
separation through the parameter Ag, and through the a-dependent quantity r in the lower 
integration limits in Eqs. (ITT)) and (13"7|) . As for the correction factor k defined in Eq. (1401) . 
it depends on separation only through the parameter Ag. 

The computational results for ft( nm ) (lines 1-6) and k (line 7) at T = 300 K are presented 
in Table I. Column 1 indicates separations a at which the values of 

K (nm) are 

computed (from 
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100 to 500 nm). Columns 2, 3, 4, and 5 contain the results obtained for Aq = 0.01, 0.05, 0.1, 
and 0.5, respectively. Keeping in mind the definition of Aq, the value of 9 = 2aAg/L in each 
column changes with the change of separation. To determine the application region of the 
multiplicative approach, we assume that this approach should not lead to an error larger 
than the error due to the PFA. For the parameters of a cylinder we take L = R = 100 /iin. 
Then, for instance, at a = 100 nm the relative error of the PFA is equal to 0.3a/ R = 0.0003 
(see Sec. II). If A e = 0.01 (see the column 2 of Table I) we have k - (100 nm) = 0.00006, 
i.e., much less than the error of the PFA. Keeping in mind that at a = 100 nm the inequality 
Aq < 0.01 is certainly valid at angles 6 < 4", one arrives at the conclusion that with 
these angles between a cylinder and a plate the multiplicative approach is applicable. It 
is very hard, however, to ensure angles 6 < 4" w 20 /zrad experimentally. Because of 
this, at short separations the nonmultiplicative approach should be used for the comparison 
between the experimental data and theory. With Aq = 0.01 the multiplicative approach 
is applicable at larger separations, as well as at even larger angles. As one more example 
we consider Aq = 0.05 (the column 3 of Table I). Here, k — /t^ nm ^(100nm) = 0.0015, i.e., 
larger than the error of the PFA at a = 100 nm. Thus, the multiplicative approach is 
not applicable. However, at a = 500 nm the error of the PFA is equal to 0.0015, whereas 
k — ft( nm ) (500 nm) = 0.0006. Thus, at this separation with chosen parameters of a cylinder 
the multiplicative approach is applicable at 9 < 1.7'. At larger angles the nonmultiplicative 
results in Eqs. ( )37l) and ( 1421) should be used to compute the Casimir force and its gradient 
with sufficient precision. 

V. ERROR DUE TO FINITE LENGTH OF A CYLINDER 

The Lifshitz-type formulas ( 1TT1) and ( 1TT1) were obtained under an assumption that the 
length L of the cylinder is infinitely large. In fact these formulas define the Casimir force 
and its gradient per unit length of a cylinder. Keeping in mind that in our experimental 
setup the length of the cylinder L is not only finite but might be of the same order as the 
cylinder radius R, it is important to estimate the possible error due to the finiteness of L. 
Here, we solve this problem for an ideal metal cylinder above an ideal metal plate. The 
derived results can serve as a reasonable estimation for real metal bodies as well. 

We start with the case when the plate is sufficiently large that the projection of the 



18 



cylinder on the plane (x,y) is situated within its interior [see Fig. 3(a)]. This is the case for 
the proposed experiment (see Sec. VI). It is assumed that the length of plane plates with an 
area Ldx, which replace a cylinder when we apply the PFA to calculate the Casimir force 
(see Sec. II), is finite. As a result, two real edges of length dx each are formed. Because 
of this, it is not legitimate to use either the original Casimir expression for the force per 
unit area — it 2 he /(240a 4 ), or its Lifshitz counterpart ([3]) applicable to real metals, in order 



to perform calculations in the framework of the PFA. Instead, we apply the expression [52] 
for the Casimir force acting between an infinite plate and a finite plate parallel to it of area 
S with edge length I spaced z apart, obtained using the method of world line numerics 

i7 i n\ ir^HcS hcl 

Frf( - 0) = -480?- (44) 

This equation with 7 a = 5.23 x was obtained for a scalar field with Dirichlet boundary 
conditions on the plates at zero temperature. Keeping in mind applications to the electro- 
magnetic field, the Casimir force acting between the partial plate of area Ldx, replacing the 
cylindrical surface in the PFA, and the plane plate is given by 

dFlf(z, 0) = —^T ~ 4 7a -r , (45) 

where z = z(x) is defined in Eq. ([[]), and the length of real edges is I = ldx. 

The integration of Eq. (1451) over the lower half of a cylinder in the leading order of the 
small parameter a/R results in 

C(M)= f R dF™(z,0) (46) 

J-R 



7i d hcL / R d,Ti^ a hc 



384a 3 V 2a a? 
F IM (a,0) (l + Cedf) 




where F (a, 0) is defined in Eq. (1131) and 

C ed = « 0.610. (47) 

In a similar way the gradient of the Casimir force with account of edge effects is given by 



da da V L 
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where <9F IM (a, 0)/da is defined in Eq. (!T8|) and 

C ed = ^C cd « 0.436. (49) 

We are now in a position to estimate the total relative error in the Casimir force and 
its gradient, originating from the application of the PFA to a cylinder of finite length. For 
the Casimir force this is a combination of the systematic errors due to the application of 
the PFA to an infinite cylinder, —C ex a/R, estimated in Sec. II, and that due to the edge 
effects, C e d(i/L. Taking the most conservative approach, we consider these systematic errors 
as random quantities characterized by a uniform distribution. Then the total error of the 
Casimir force determined at a 95% confidence level can be obtained from the following 
statistical ralefl 



a „ a I „„ a 2 „ n a? 



6? = min I C«- + C cdI , l.iy C e 2 x — + Cl AT2 J . (50) 

Taking R = L = 100 /im we find from Eq. §0D that 8$* « 0.07% at a = 100 nm and 
o> ot w 0.37% at a = 500 nm. 

To estimate the total error of the gradient of the Casimir force between a cylinder and a 
plate, one should use the exact expression following from Eq. flTJ 

da ~ 768a 4 V 2a V ex L, 
Cex = 5c« « 0.2062. (51) 

This shows that the error in the application of the PFA to the gradient of the Casimir force 
between an infinite cylinder and a plate is less than in the application of this method to the 
force. Combining the errors C cx a/R from Eq. ( 15T|) and C e da/L defined in Eq. ( 1491) with the 
help of the same statistical rule fl50|) . one obtains S^p, da ~ 0.05% and 0.26% at a = 100 nm 
and 500 nm, respectively. 

In the dynamic experiments measuring the gradient of the Casimir force between a sphere 
and a plate by means of a micromachined oscillator 
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36| the Au-coated sphere was at- 



tached to the optical fiber. As a result, the projection of the part of the sphere near the 
point of attachment on the plane of an oscillator plate was situated beyond the plate edge. 
As was noted above, in the proposed experiment the cylinder is attached to an oscillator 
and the cylinder's projection on the plate belongs to its interior (see Sec. VI). However, in 
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some other experimental configurations it may happen that the projection of the part of the 
cylinder is beyond the edge of the plate. We now consider the role of edge effects in the case 
when the projection of a cylinder axis on the plate is at distance L\ < R from the edge of a 
plate [see Fig. 3(b)]. In this case one half of the cylinder is completely above the plate. As 
to the other half, only the portion x < L\ is above the plate. 

We next obtain the expression for the Casimir force taking into account edge effects, 



where dF^ (z, 0) is defined in Eq. <$E§, and H = R - ^R 2 - L\ [see Fig. 3(b)]. In the 
following we assume L%, H 3> a. Here, the upper integration limit L\ [instead of R in 
Eq. f|46|) ] takes into account that only part of the cylinder is above the plate. The second 
term on the right-hand side of Eq. (152]) is due to the presence of an additional real boundary 
of length L. This boundary belongs to one of the partial plane plates, with area Ldx 
replacing a cylinder when the PEA is used which is situated at a height a + H above the 
edge of an oscillator plate. Equation ( 1521) can be equivalently presented in the form 



Performing the integrations and retaining only the leading terms in the small parameters 
a/R and a/H we arrive at the result 




(52) 




(53) 



F e I d M (a,0)=F IM (a,0) 1 + C ed 



a 



(54) 



L 




Here, the following notations are introduced: 




(55) 



where 



fx{z) = 120 - 160z 2 + 35^ 4 + 13z 



(56) 



+ 4Vl-z 2 (30 - 27 z 2 -z 4 - 2z 6 ), 



f 2 ( z ) = 4 - 5z 2 + z 5 + VT^(4: - 3z 2 - z A ). 
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If L\ = R (i.e., the cylinder is completely over the plate) we have /i(l) = ^2(1) = /(l) = 0, 
and Eq. ( I54p reduces to Eq. (j4U|) . [Note that in this case the last term on the right-hand 
side of Eq. (I54"|) is of order a 3 /R 3 and, thus, should be disregarded in the framework of the 
PFA.] 

It is easily seen that two additional terms on the right-hand side of Eq. ( I54p are negligibly 
small. For example, for the same parameters R = L = 100 //m, as above, and L\ = 25 /im the 
contribution of these terms does not exceed 10~ 4 % at a = 100 nm and 0.05% at a = 500 nm. 
With the increase of L\ the contribution of the additional terms in Eq. (154")) further decreases. 
For instance, for Li = 50 /mi it does not exceed 5 x 10 _7 % at a = 100 nm and 1.5 x 10 _4 % 
at a = 500 nm. 



VI. PROPOSED EXPERIMENTAL SETUP 

In this section we propose an experiment to measure the Casimir interaction between 
a cylinder and a plate using an approach similar to what was used in previous studies for 



sphere-plate geometry j32|-|36|, |70j . A sensitive micromechanical torsional oscillator (MTO) 
is at the heart of the setup since it allows detection of the Casimir force in a very precise 
manner. An incorporated optical fiber two-color interferometer allows the measurement of 
the fiber-to-substrate separation at all times. As it was done in our previous experiments, 
an electrostatic calibration will allow us to determine the interaction force, as well as the 
separation a. 

We propose to integrate a metallic cylinder onto a MTO and to detect the deflections 
experienced by the micromechanical device when a metallic plate approaches the cylinder 
from above (see Fig. 4 for a schematic of the setup). Section II presents theoretical results 
for the Casimir force and for the gradient of the Casimir force between a cylinder and a 
plate. Both these quantities can be measured using the suggested setup. The Casimir 
force is measured in the static mode 32]. The gradient of the Casimir force is measured in 



the dynamic mode, where the separation distance between a cylinder and a plate is varied 



harmonically at the resonant frequency of the MTO [33j . The experimental setup allows for 
a force sensitivity of ~ 0.1 pN and a determination of the resonant frequency of the MTO 
to better than 6mHz. In the dynamic mode, the gradient of the Casimir force is obtained 
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from the resonant frequency of the oscillator 

Wres = H T^J- (57) 

Here, ujq is the natural angular resonance frequency of the oscillator, co> res is the angular 
resonant frequency of the oscillator in the presence of the Casimir force, b is the lever arm 
of the interaction, and I is the moment of inertia of the oscillator. 

In order to minimize errors associated with misalignments between the torsional axis 
of the MTO and the cylinder's axis, we propose to integrate the metallic cylinder onto 
the MTO using state-of-the-art semiconducting processing. Currently available lithographic 
steppers allow alignment between different structural layers in a mechanical device with an 
accuracy of 20 nm or better. By using a monolithic fabrication process it would be possible 
to directly fabricate a metallic cylinder on top of a MTO similar to the one previously used. 
In order to fabricate a cylinder with a pre-determined cylindrical shape and smooth surface 
it will be necessary to use grayscale masks or laser based maskless lithography technology 



71j . These technologies permit the creation of three-dimensional structures with smooth 



surfaces and precise critical dimensions. A simpler method to fabricate a cylindrical object 
on top of a MTO involves deposition and patterning of photosensitive polymers (photoresist), 
and subsequent thermal reflow of the polymer into a curved shape. This technique relies on 
surface tension and plastic flow of photoresist and is used to produce smooth curved objects. 
Reflow techniques are much simpler to implement than grayscale lithography techniques, 
but they are much more limited regarding the final shape of the cylindrical objects. The 
photoresist shape can be transferred to the MTO using conventional dry etching techniques 
like reactive ion etching (RIE). Once the cylindrical shape has been transferred to the MTO 
it can be coated with a metallic film by evaporation or sputtering methods. 

The actual cylindrical shape z = f(x, y) can be measured very precisely, better than 1 nm 
error in z, by using a non-contact optical profilometer. Once this shape is known, it can be 
introduced in Eq. ([I]) to calculate the Casimir force acting on the cylinder, as described in 
Sec. II. 

During the measurements, the cylinder-MTO structure will be brought close to the Au- 
coated plane using a 5-axis manipulator with an angular error of 10~ 5 radians. This angular 
error allows for an alignment of the cylinder and the plate to within the ranges calculated 
in Sec. IV. 
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Further improvements in the alignment between the axis of the MTO and the plate can 
be obtained by sacrificing the ability to perform continuous measurements as a function 
of separation. In this approach the oscillator/cylinder unit is manufactured as previously 
described, but the Au-coated plate on a Si substrate is engineered with three equal pillars of 
known height. The MTO chip is positioned on top of the pillars and a measurement of the 
Casimir interaction at a fixed separation can be performed. By changing the height of the 
pillars it is feasible to reproduce a curve of the Casimir interaction as a function of separation. 
Since the height of the pillars can be controlled to better than 10 nm, misalignments smaller 
than 1 /zrad are possible. 



VII. CONCLUSIONS AND DISCUSSION 



In the foregoing we have investigated the possibilities of measuring the thermal Casimir 
force and its gradient between a plate and a microfabricated cylinder attached to a microma- 
chined oscillator. Using the PFA we have derived the Lifshitz-type formulas for the thermal 
Casimir force for a cylinder and a plate made of real materials. From the comparison with 
exact results available for an ideal metal cylinder above a plate and for both ideal and real 
metals in the configuration of a sphere above a plate, we have estimated the error resulting 
from the use of the PFA. It turns out that for reasonable experimental parameters this error 
does not exceed a small fraction of a percent. In the limiting case of high temperature (large 
separation) the analytic expressions for the Casimir force and its gradient are obtained. 

The derived Lifshitz-type formulas were used to perform numerical computations of the 
thermal correction to both the Casimir force and its gradient. In so doing the two alternative 
theoretical approaches proposed in the literature (the Drude model and the plasma model 
approaches) have been used. The predicted results of both static and dynamic experiments 
for the measurement of the thermal correction were discussed. Specifically, it was concluded 
that the thermal correction, as predicted by the Drude model approach, is large enough 
to be measurable at short separations below 1 /im. With respect to the thermal correction 
predicted by the plasma model approach, it can be measured only at separations of a few 
micrometers. The performed computations, when compared with the experimental data, 
may shed additional light on the experimental exclusion of the Drude model approach in 
Refs. 132 



361 ] which was made using only the configuration of a sphere above a plate. 
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Keeping in mind the applications of the developed formalism for comparison with exper- 
imental data, we have derived expressions for the thermal Casimir force and its gradient 
taking into account the nonparallelity of the cylinder and plate. This was done in a non- 
multiplicative manner including the role of correlations between geometrical and material 
properties. In the specific case of ideal metal bodies at zero temperature, the previously 
obtained results were confirmed. We have also performed numerical computations of a 
correction to the Casimir force due to nonparallelity of a cylinder and a plate using both 
nonmultiplicative and multiplicative approaches, and determined the application region of 
the latter. 

Taking into account that microfabricated cylinders to be used in experiments are re- 
stricted in length (which may be of the order of the cylinder radius), we have estimated 
errors in computations arising from the finiteness of a cylinder. The two experimentally 
relevant cases were considered when the cylinder is completely over the plate or is partially 
beyond the plate edge. In both cases it was shown that for typical experimental parameters 
the total error due to application of the PEA and due to edge effects is sufficiently small, 
and does not depend on the ratio of the cylinder radius to the cylinder length. This opens 
opportunities for using cylinders with L m R, with no loss in the accuracy of the theoretical 
expressions derived for infinitely long cylinders. 

We have presented a scheme of the proposed experiment for measuring thermal Casimir 
force and its gradient in the configuration of an Au-coated plate and a microfabricated 
cylinder attached to the micromachined oscillator. Both the static and dynamic versions of 
this experiment are considered. The values of main experimental parameters that are aimed 
to achieve are estimated. The proposed experiment is promising as a source of additional 
information on the problem of thermal Casimir force. 
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FIG. 1: (Color online) The relative thermal corrections to the Casimir force and its gradient in 
cylinder-plate configuration computed using the Drude model approach at T = 300 K (lines 1 and 
2, respectively) are shown as functions of separation in the region from 0.1 to 5//m. 
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FIG. 2: (Color online) The relative thermal corrections to the Casimir force and its gradient in 
cylinder-plate configuration computed using the plasma model approach at T = 300 K (lines 1 and 
2, respectively) are shown as functions of separation in the region (a) from 0.1 to 5 /xm and (b) 
from 0.1 to 1 /im. 
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FIG. 3: The cylinder and the plate in cases when (a) the projection of the cylinder on the plane 
(x, y) belongs to the interior of the plate and (b) part of the projection of the cylinder on the plane 
(x, y) is beyond the edge of the plate. 
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FIG. 4: (Color online) Schematic of the experimental setup (see text for further discussion). The 
figure is not to scale. 
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TABLE I: The lines from 1 to 6 contain the values of the quantity k*" 111 ) quantifying the effect of 
the nonzero angle between a cylinder and a plate on the Casimir force within the nonmultiplicative 
approach. In line 7 the correction factor n taking the nonzero angle in a multiplicative way is 
presented. All results are computed at different values of the parameter Ag = 9L/(2a). 







Ag = 0.01 


Ag = 0.05 


Ag = 0.1 


Ag = 0.5 


K (nm < 


(100 nm, Ag) 


1.00020 


1.0051 


1.0207 


1.7888 


K (nm] 


(150 nm, 4*) 


1.00021 


1.0053 


1.0215 


1.8230 


K (nm] 


(200 nm,^) 


1.00022 


1.0055 


1.0221 


1.8526 


K (nm] 


(300 nm, Ag) 


1.00023 


1.0057 


1.0231 


1.8997 




(400 nm, Ag) 


1.00023 


1.0058 


1.0237 


1.9335 




(500 nm, Ag) 


1.00024 


1.0060 


1.0242 


1.9585 




K{Ag) 


1.00026 


1.0066 


1.0267 


2.1176 
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